New tests for a singularity of ideal MHD 
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Analysis using new calculations with 3 times the resolu- 
tion of the earlier linked magnetic flux tubes confirms the 
transition from singular to saturated growth rate reported by 
Grauer and Marliani || for the incompressible cases is con- 
firmed. However, all of the secondary tests point to a transi- 
tion back to stronger growth rate at a different location at late 
times. Similar problems in ideal hydrodynamics are discussed, 
pointing out that initial negative results eventually led to bet- 
ter initial conditions that did show evidence for a singularity 
of Euler. Whether singular or near-singular growth in ideal 
MHD is eventually shown, this study could have bearing on 
fast magnetic reconnection, high energy particle production 
and coronal heating. 



The issue currently leading to conflicting conclusions 
about ideal 3D, incompressible MHD is similar |QJ|] to 
what led to conflicting results on whether there is a sin- 
gularity of the 3D incompressible Euler. With numeri- 
cal simulations, it was first concluded that uniform mesh 
calculations with symmetric initial conditions such as 3D 
Taylor- Green were not yet singular Next, a prelimi- 
nary spectral calculation Q] found weak evidence in favor 
a singularity in a series of Navier-Stokes simulations at 
increasing Reynolds numbers, but larger adaptive mesh 
or refined mesh calculations did not support this result 
|^],[|. Eventually, numerical evidence in favor of a sin- 
gularity of Euler was obtained using several independent 
tests applied to highly resolved, refined mesh calculations 
of the evolution of two anti-parallel vortex tubes . To 
date, these calculations have met every analytic test for 
whether there could be a singularity of Euler. 

Several other calculations have also claimed numerical 
evidence for a sing ularity of Euler While in all of 

these cases the evidence is plausible, with the perturbed 
cylindrical shear flow JT3] using theBKM Halloo test @, 
for none has the entire battery of tests used for the anti- 
parallel case been applied. We have recently repeated 
one of the orthogonal cases S and have applied the BKM 
test successfully. In all cases using the BKM test, \uj\\oo ~ 
A/(T C - t) with A w 19. 

To be able to make a convincing case for the existence 
of a singularity in higher dimensional partial differential 
equations, great care must be taken with initial condi- 
tions, demonstrating numerical convergence, and com- 
parisons to all known analytic or empirical tests. On the 
other hand, if no singularity is suspected, some quan- 
tity that clearly saturates should be demonstrated, such 
as the strain causing vorticity growth ||. It is an even 
more delicate matter to claim that someone else's calcu- 



lations or conclusions are incorrect. If it is a matter of 
suspecting there is inadequate resolution, one must at- 
tempt to reproduce the suspicious calculations as nearly 
as possible and show where inadequate resolution begins 
to corrupt the calculations and how improved resolution 
changes the results. 

An example of how a detailed search for numerical 
errors should be conducted can be found in the exten- 
sive conference proceeding ]T^ ] that appeared prior to 
the publication of the major results supporting the ex- 
istence of a singularity of Euler for anti-parallel vortex 
tubes 0. The primary difference with earlier work was 
in the initial conditions. It was found that compact pro- 
files fl3f| were an improvement, but only if used in con- 
junction with a high wavenumber filter. Otherwise, the 
initial unfiltered energy spectrum of the bent anti-parallel 
vortex tubes went as k~ 2 . Oscillations in the spectrum 
at high wavenumber in unfiltered initial conditions for 
linked magnetic flux tubes are shown in Figure [I], show- 
ing that the initial MHD spectrum is steep enough that 
eventually these oscillations are not important. 

Filtered and unfiltered 
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FIG. 1. Filtered and unfiltered initial and final spectrum . 
The unfiltered spectrum is initialized on a 384 3 mesh. 

The purpose of this letter is to address the claim that 
a new adaptive mesh refinement (AMR) calculation by 
Grauer and Marliani Q supercedes our uniform mesh cal- 
culations 0| and that eventually there is a transition to 
exponential growth. Note that this claim was made with- 
out any evidence for whether their numerical method was 
converged. In all of our earlier calculations, once the cal- 
culations become underresolved, we also saw transitions 
to exponential growth. 

Not knowing exactly the initial condition used by the 
new AMR calculations (|] , where and how much grid re- 
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finement was used, and the short notice we have been 
given to reply has proven a challenge. Fortunately, we 
were in the process of new 648 3 calculations in a smaller 
domain of 4.3 3 , yielding effectively 3 times the local res- 
olution of our earlier work Jj| in a (2ir) 3 domain on a 
384 3 mesh. The case with an initial flux tube diameter 
of d = 0.65, so that the tubes slightly overlap, appears to 
be closer to their initial condition and so will be the focus 
of this letter. The importance of our other initial condi- 
tion, with d = 0.5, and no initial overlap of the tubes, is 
that it is less influenced by an initial current sheet that 
forms near the origin and is claimed to be the source of 
the saturation of the nonlinear terms. This was used for 
the compressible calculations. 



shows that either an exponential or a singular 1/(T C — t) 
form could fit the data, while the inset shows that taking 
an estimated singular time of T c = 2.15 and multiplying 
by (T c — T) that at least || J||oc and |M|oo have consistent 
singular behavior over this time span. The strong growth 
of Pnj = J dVfaeijWj - LOidi 3 Jj - 2s i3k Jidjeetk), which 
is the production of J dV(ui 2 + J 2 ), is discussed below. 
The 384 3 curve for 1/|| J| demonstrates that lack of res- 
olution tends to exaggerate exponential growth. For the 
compressible calculations it can be seen that there also 
is an exponential regime that changes into a regime with 

l/||J||oo~(T c -t). 
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FIG. 2. Replot of ||J|| oo ; W CO) and Pnj for the new in- 
compressible calculations on 4.3 3 domain with initial condi- 
tion d — 0.65 in semi-log coordinates. All plots are from the 
648 3 calculation except one 384 3 plot of 1/||J||. Exponential 
and inverse linear fits are shown for t = 1.75 to 1.98. Each 
works equally well for || J||oo, inverse linear is better for ||w||aoj 
and exponential is better for Pnj- Multiplying ||J||oo, IMloo, 
and Pnj by (T c — t) in the inset emphasizes that || J||oo and 
M|oo might be showing consistent singular behavior. The 
large growth of Pnj is discussed. 

Using semi-log coordinates, Figure || plots the growth 
of H^Hoo and 1 1 J | |oo f° r our new high resolution incom- 
pressible calculation and Figure || plots ||J||oo for a new 
compressible calculation. By taking the last time all rel- 
evant quantities on the 384 3 and 648 3 grids were con- 
verged, || J|| being the worst, then by assuming that the 
smallest scales are decreasing linearly towards a possible 
singular time, an estimate of the last time the 648 3 cal- 
culation was valid was made. To test exponential versus 
inverse linear growth, fits were taken between T = 1.72 
and 1.87, then extrapolated to large T. The large figure 
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FIG. 3. Semi-logarithmic plot of ||J||oo for a compressible 
240 3 calculation in a domain of size 4 (dotted line: filtered, 
and solid line: unfiltered initial conditions) together with fits 
to exponential growth and blow-up behavior, respectively. 
The latter are better fits at later times. 



Using the new incompressible calculations and apply- 
ing the entire battery of tests, based upon Figure E we 
would agree that for the incompressible case there is a 
transition as reported j|] and signs of saturation at this 
stage are shown below. Whether the transition is to ex- 
ponential for all times as claimed [|| , or whether there is 
a still later transition to different singular behavior, will 
be the focus of this letter. We will look more closely at 
the structure of the current sheet we all agree exists 
for signs of saturation. 

The case against a singularity in early calculations of 
Euler [p| p^Jl5[ | was the appearance of vortex sheets, and 
through analogies with the current in 2D ideal MHD, a 
suggestion that this leads to a depletion of nonlinearity. 
The fluid flow most relevant to the linked flux rings is 
3D Taylor-Green, due to the initial symmetries jjj. For 
both TG and linked flux tubes, two sets of anti-parallel 
vortex pairs form that are skewed with respect to each 
other and are colliding. In TG, just after the anti-parallel 
vortex tubes form there is a period of slightly singular 
development. This is suppressed once the pairs collide 
with each other, and then vortex sheets dominate for a 
period. The vortex sheets are very thin, but go across 
the domain, so fine localized resolution might not be an 
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advantage at this stage. At late phases in TG, the ends 
of the colliding pairs begin to interact with each other, 
so that at 4 corners locally orthogonal vortices begin to 
form. Due to resolution limitations, an Euler calculation 
of Taylor-Green has not been continued long enough to 
determine whether, during this phase, singular behavior 
might develop. We would draw a similar conclusion for 
all of MHD cases studied to date P, |l6| , pT| , that there 
might not be enough local resolution to draw any final 
conclusions even if AMR is applied. 

While Taylor-Green has not been continued far enough 
to rule out singularities, the final arrangement of vortex 
structures led first to studies of interacting orthogonal 
vortices || , and then anti-parallel vortices (see references 
in ) . Both of these initial conditions now appear to de- 
velop singular behavior. An important piece of evidence 
for a singularity of Euler was that near the point of a 
possible singularity, the structure could not be described 
simply as a vortex sheet. Therefore, there is a precedent 
to earlier work suggesting sheets, suppression of nonlin- 
earity, and no singularities to later work showing fully 
three-dimensional structure and singular behavior. 

The initial singular growth of ||J||oo and |M|oo for 
the linked flux rings, then the transition to a saturated 
growth rate, might be due to the same skewed, anti- 
parallel vortex pair interaction as in Taylor-Green. Even 
if this is all that is happening, the strong initial vortic- 
ity production and shorter dynamical timescale (order of 
a few Alfven times) than earlier magnetic reconnection 
simulations with anti-parallel flux tubes |l7f is a signif- 
icant success of these simulations. It might be that the 
vortices that have been generated are strong enough to 
develop their own Euler singularity. However, the in- 
teresting physics is how the magnetic field and current 
interact with the vorticity. Do they suppress the ten- 
dency of the vorticity to become singular, or augment 
that tendency? 
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FIG. 4. Positions of || and \\u\loo for d = 0.65 in a 4.3 3 
domain. 

One sign for saturation of the linked flux ring inter- 



action would be if the strongest current remains at the 
origin in this sheet. Figure f| plots the positions of | J||oo 
and ||w||oo from the origin as a function of time. Dur- 
ing the period where exponential growth is claimed [^), 
|| J | |oo is at the origin, which would support the claims of 
saturation. However, this situation does not persist. 

By analogy to the movement of the L°° norms of the 
components of the stress tensor itjj in Euler, we expect 
that the positions of \\J\\oo and ||w||oo should approach 
each other and an extrapolated singular point in ideal 
MHD. Figure | supports the prediction that the positions 
of 1 1 J | |oo an d IMIoc should approach each other but so far 
not in a convincingly linear fashion. This is addressed 
next. We have similar trends for the positions of || J||oo 
and 1 1 oj ||oo in the compressible calculations. 




FIG. 5. For t — 1.97 on the inner 162 3 grid points, the cur- 
rent sheet is shown with arrows of Tf overlaid in dark. The 
current through the (x/y = z) plane containing ||J||oo is in 
lower right. Contours of |J| 4 are shown to emphasize where 
||J||oo is located. Dark lines are Tl and light lines are ~uf 
that originated in the vicinity of || J||oo- The vortex lines are 
predominantly those in the double vortex rings that were orig- 
inally generated by the Lorenz force, then became responsible 
for spreading out the current sheet. Where the lines cross 
in the upper left and lower right corners are around the lo- 
cations of ||J||oo, which due to symmetries are different views 
of the same structure. Near ||J||oo, B nearly overlies and is 
parallel to ~ut and both and ~uf are nearly orthogonal to 
their partners across the current sheet, where and ~ut are 
anti-parallel. Taken from the d = 0.65 calculation in a 4.3 3 
domain on a 648 3 mesh. 

Figure gives an overall view of the current, vorticity 
and magnetic field around the inner current sheet. The 
vortex pattern has developed out of the four initial vor- 
tices, two sets of orthogonal, anti-parallel pairs that are 
responsible for the initial compression and stretching of 
the current sheet. By this late time, the ends of the those 
vortices have begun to interact as new sets of orthogonal 
vortex pairs. The lower right inset in Figure [| is a 2D 
x/(y = z) slice through this domain that goes through 
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||«7||oo at t — 1-97 to show that while || J||oo is large at the 
origin (0,0,0), || J\\oo is larger where it is being squeezed 
between the new orthogonal vortices. Along one of the 
new vortices B is parallel to and overlying a? and on the 
orthogonal partner they are anti-parallel and overlying. 

The location of ||w||oo is not in the vortex lines shown, 
but is on the outer edges of the current sheet. Therefore, 
the exact position of |M|oo in Figure [| is an artifact of 
the initial development and does not accurately reflect 
the position of ~uf most directly involved in amplifying 
|| J | |oo, which is probably why the positions of \\J\\oo and 
H^Hoo are not approaching each other faster. The con- 
tinuing effects of the initial current sheet is probably also 
behind the strong exponential growth of Pjij in Figure ||, 
stronger even than the the possible singular growth of 
|| J | |oo and |M|oo in the inset. More detailed analysis in 
progress should show that near the position of ||</||oo, the 
growth of Pqj and the position of ||w||oo are more con- 
sistent with our expectations for singular growth and has 
already shown that some of the components of Pjij have 
consistent singular growth. 

As noted, for Euler all available calculations find 
Mloo ~ A/(T C — t) with A a 19. A represents how 
much smaller the strain along ||w||oo is than H^Hoo. Here, 
A 4, indicating stronger growth in ||w||oo for ideal 
MHD than Euler. Another Euler result was that the 
asymptotic energy spectrum as the possible singularity 
was approached was fc~ 3 , whereas purely sheet-like struc- 
tures in vorticity should yield fc~ 4 spectrum. k~ 3 indi- 
cates a more complicated 3D structure than sheets. In 
Figure [I] the late time spectra are again k~ 3 . 

The next the initial condition we will investigate will be 
magnetic flux and vortex tubes that nearly overlay each 
other and are orthogonal to their partners. Our new cal- 
culations of orthogonal vortex tubes for Euler show that 
they start becoming singular as filaments are pulled off 
of the original tubes and these filaments become anti- 
parallel, suggesting that the fundamental singular inter- 
action in Euler is between anti-parallel vortices. Whether 
the next step for ideal MHD is to become anti-parallel 
or something else can only be determined by new cal- 
culations. AMR might be useful, but great care must 
be taken with the placement of the inner domains and a 
large mesh will still be necessary. The complicated struc- 
tures in the domain in Figure || are not fully contained 
in this innermost 162 3 mesh points and the innermost 
domain should go out to the order of 300 3 points. There 
are examples of how to use AMR when there are strong 
shears on the boundaries of sharp structures This 
uncertainty of where to place the mesh is why we believe 
in using uniform mesh calculations as an unbiased first 
look at the problem. 

These final results are hardly robust and their useful- 
ness is primarily to suggest a new more localized initial 
condition and to show that none of the calculations to 
date is the full story. For J and u> to show singular be- 
havior as long as they have has been surprising. Recall 
that for Euler, velocity, vorticity and strain are all man- 



ifestations of the same vector field, but for ideal MHD 
there are two independent vector fields even though the 
only analytic result in 3D is a condition on the combi- 
nation, /dV[|M| 00 (t) + ||J||oo(i)]d£ -> oo (19). Even- 
tually, one piece of evidence for singular growth must 
be a demonstration of strong coupling between the cur- 
rent and vorticity so that they are acting as one vector 
field. It could be that our strong growth is due to the 
strongly helical initial conditions and there are no singu- 
larities. This would still be physically interesting since 
helical conditions could be set up by footpoint motion in 
the corona. 

Could the magnetic and electric fields blow up too? 
There are signs this might be developing around the fi- 
nal position of || J||oo, in which case there might exist a 
mechanism for the direct acceleration of high energy par- 
ticles. This has been considered on larger scale |2fJ, but 
to our knowledge a mechanism for small-scale production 
of super-Dreicer electric fields has not been proposed be- 
fore. A singular rise in electric fields could explain the 
sharp rise times in X-ray production in solar coronal mea- 
surements [ pl[ , which could be a consequence of particle 
acceleration coming from reconnection. This would also 
have implications for the heating of the solar corona by 
nanoflares ]22| and the production of cosmic rays. 
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National Science Foundation. 
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